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Nucleation on top of two-dimensional islands with step edge barriers is investigated using scaling 
arguments. The nucleation rate is expressed in terms of three basic time scales: The time interval 
between deposition events, the residence time of atoms on the island, and the encounter time required 
for i* + 1 atoms forming a stable nucleus to meet. Application to the problem of second-layer 
nucleation on growing first layer islands yields a sequence of scaling regimes with different scaling 
exponents relating the critical island size, at which nucleation takes place, to the diffusion and 
deposition rates. Second layer nucleation is fluctuation-dominated, in the sense that the typical 
number of atoms on the island is small compared to i* -I- f , when the first layer island density 
exponent x satisfies -I- f ) < 2. The upper critical nucleus size, above which the conventional 
mean-field theory of second layer nucleation is valid, increases with decreasing dimensionality. In 
the related case of nucleation on top of multilayer mounds fluctuation-dominated and mean-field 
like regimes coexist for arbitrary values of the critical nucleus size i* . 
PACS: 68.35.Fx, 81.15.Aa, 05.40.-a 



I. INTRODUCTION 



Atomistic nucleation theory was developed in the sixties and seventies in response to the realization that critical 
nuclei of atomic dimensions are common in far from equilibrium thin film growth During the past decade the 

subject has experienced a vigorous renaissance, driven by the invention of atomic scale imaging techniques as well as 
by the availability of realistic large scale simulations Q . While this has resulted in a wide range of new theoretical 
developments, the classical theory has stood the test of time remarkably well. This is surprising because classical 
nucleation theory relies, in the jargon of statistical mechanics, on an approximation of mean field type, in the sense 
that it considers only the evolution of spatially averaged concentrations of the different surface species - adatoms and 
clusters ~ and thereby neglects effects of spatial correlations and fluctuations. 

It was recently pointed out that the classical theory manifestly fails in describing nucleation events on top of islands 
bounded by strong step edge barriers |^J^. The onset of second layer nucleation on top of the first layer islands 
is an important factor determining the growth mode and morphology of multilayer films. Tersoff and coworkers 0| 
developed a theory of second layer nucleation based on the classical mean field estimate 

a;,nf ~ i^An''+^ (1) 

for the rate uj of nucleation events on the island, in terms of the in-layer hopping rate i', the island area A, the adatom 
density n and the size i* of the critical (= largest unstable) cluster. In Refs. it was shown that for «* = 1, 2 
the mean field estimate disagrees with microscopic considerations ]^-[To[| and computer simulations. The failure of 
mean field theory was qualitatively attributed to the small number of atoms typically present on the island, and the 
associated large fluctuations in the occupancy of the island. 

The question then arises whether mean field behavior may be recovered at larger values of the critical nucleus 
size, where nucleation necessarily involves a greater number of atoms, and thus fluctuation effects should be reduced. 
A partial answer to this question was provided in recent work by Heinrichs, Rottler and Maass (HRM) |ll[]. They 
distinguish between two different nucleation modes: In the fluctuation- dominated mode the mean number of adatoms 
on the island (in the quasi-stationary state prior to nucleation) is much smaller than the number i* + 1 required for 
nucleation, and therefore nucleation is a rare, large fluctuation, while in the mean field regime, which corresponds to 
the classical approach of 0, the mean number of adatoms exceeds i* + 1. Consistency considerations then show that 
second layer nucleation on top of first layer islands is fluctuation-dominated for i* < 2 and mean-field-likc for i* > 3. 

In this paper I elaborate on the observations of HRM, with the goal of further clarifying the role of fluctuations 
in second layer nucleation. I begin by rederiving the main results of HRM within the approach of 1^, which relies 
on expressing the nucleation rate in terms of the basic time scales of the process. While less quantitative than the 
treatment of this provides, perhaps, a somewhat more transparent and unified view. A simple, universal criterion 
for fluctuation-dominated second layer nucleation is found which involves the first layer island density exponent and 
the critical nucleus size (Eq.(p8|)). I then extend the analysis to d-dimensional nucleation, finding the same general 
structure as in the two-dimensional case, but with a decreasing influence of fluctuations with increasing d; as a 
byproduct, I recover the expression for the one-dimensional first layer island density exponent recently derived in Rcf. 
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|_2|. Finally, I generalize the theory of nucleation on top of multilayer mounds developed in |^ to i* > 1. In contrast 
to the case of nucleation on first layer islands, here fluctuation-dominated and mean field like nucleation regimes are 
found to coexist for arbitrary values of i*. This illustrates that the importance of fluctuations depends crucially on 
the growth situation. 



II. STATIONARY RATES OF FLUCTUATION-DOMINATED NUCLEATION 



We consider a two-dimensional island of fixed linear extension L. Atoms are deposited onto the island at rate F, 
they migrate on the island with an in-layer hopping rate v and descend from the island at the inter-layer hopping 
rate v' . The mean time interval between subsequent arrivals of atoms on the island is then 

A^-^ (2) 
and the mean residence time of an atom on the island in steady state is of the form ]5|,|6|,pT| 

aL^ bL 



V 



(3) 



where a and h are geometry-dependent constants. The first term on the right hand side of is of the order of the 
traversal time rtr ~ I? jv required for the atom to migrate across the island, while the second term describes the 
increase in the residence time due to the step edge barriers. The second term dominates if L <C /eSj where 

Zes « vjv' (4) 

is the Ehrlich-Schwoebel length 

The stationary mean number of atoms on the island is given by r/At PJll[|, and hence the condition for the 
fluctuation-dominated mode of HRM reads 

it « 1- (5) 

It will turn out to be useful to write this in the form L <C if, where the flux scale Lp given by 



^^^^ Wf)V4 : i»/j,3. (6) 



For given i, V and v\ Eq.d) can always be satisfied by making the flux sufficiently small. 

Nucleation with a critical nucleus size i* is treated within the "noninteracting particle model" [ pl| , in which the 
lifetimes of all unstable clusters are assumed comparable to the inverse hopping rate Xjv. Then the third relevant 
time scale in the problem, in addition to and (|3|) , is given by the encounter time Tcnc required for the z* -I- 1 atoms 
to find each other on the island. In order of magnitude we have JlH 



IV 



A simple derivation is given below in Section 

Under the condition (^, each nucleation event is uniquely associated with the deposition of the i* -\- I'th atom. 
The number of nucleation events per unit time can then be written as 

W = Pnuc/At, (8) 

where pnuc is the probability for a freshly deposited atom to participate in a nucleation event. The latter is generally 
of the form 

Pnuc = K' + l ■ Pone, (9) 

where denotes the probability that i* atoms were present on the island prior to the arrival of the i* -f I'th, and 

Pcnc is the probability that the i* -I- 1 atoms encounter each other before one of them leaves the island again. The 
quantity p^.+i is given by ||] 
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The encounter probability can be estimated as p 



-r j Tcnc 



(11) 



where for simpHcity we have equated the residence time of the + 1 atoms with that of a single atom (in fact the two 
differ by a factor i* + 1). Thus penc is of order unity for r/rcnc ^ 1 and of the order penc ~ T/renc when r/renc 1. 

In the case i* = 1 the encounter time (0) is of the same order as the traversal time -^V^^i whereas the residence 
time (^ is of the order of Ttr or larger. Thus one always has Ponc — C'(l) and p„uc ~ V2 @- In contrast, for i* > 1 
it is possible to have a transition from r ^ Tcnc to r Tone with increasing island size. This obviously requires that 
J^/i^' ^ 1, so that T ^ Ljv' Tone for small islands, and occurs at the characteristic length scale 



{vjv') 



l/(2»* 



(12) 



We note that this is small compared to the Ehrlich-Schwoebel length, and therefore the transition to encounter-limited 
nucleation occurs before the step edge barriers become irrelevant, in the sense that t becomes of the same order as 

Ttr. 

In the encounter- limited regime the expression (H) agrees in order of magnitude with the mean field estimate (|^). 
To see this, note that the adatom density on the island is of order T/(AiL^) = Ft, and hence 



(L^At) 



i'+l 



At \AtJ 



Pi*+1 -Pc 

At 



(13) 



In this sense the mean field approach remains applicable also outside of its strict range of validity |11|. 

We summarize the preceding considerations in the following three expressions for the fluctuation-dominated nucle- 
ation rate, which apply with increasing island size L: 
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At \AtJ 



F{F/v 



7-3i*+2 



L ^ Ly, (regime II) 



(14) 



1 / T \ + l 



VAtJ 



Lx < i < Zes (regime III) 



(15) 



— f- 

Tone VAt 



I' + l 



F{F/iyy'L 



2i*+4 



L > Zes (regime IV). 



(16) 



Regime III does not exist for i* = 1. For i* > 1, the encounter-limited regimes III and IV extend also into the mean 
field realm, while regime II occurs only in the fluctuation-dominated nucleation mode. 

Note that Eqs.([l^) and ( p^ correspond to the regimes II and I identified in |^ for i* — 2. Here we pr efer to 
enumerate the regimes in accordance with HRM |ll[ . Their regime I will be dealt with below in Section [II B . 



III. CONSISTENT SCALING REGIMES IN TWO DIMENSIONS 



While Eqs.( |l4|Jl5| , |lq ) provide, in the specified ranges of parameters, valid expressions for the nucleation rate on top 
of an island of given size, in their application to a specific growth situation additional consistency requirements arise 
pTf : When a given expression for u is used to compute the island size at which nucleation occurs, it has to be checked 
that this island size does indeed lie in the appropriate range. 

In this section we derive the resulting scaling regimes for second layer nucleation in two dimensions. Following 
earlier work 11 1 we assume a population of equal sized first layer islands of density N. The island size increases 
with time t or coverage 6 = Ft according to the growth law 



(17) 



For a nucleation rate of the form 



uj = FQ.L'' 



(18) 



the critical island size at which a second layer nucleus forms is then given by j^] 

L^^ {Nn)-^/^''+^K (19) 

We further introduce the notation pl 

T = v/F, a = v'/v = l:^l, (20) 
where usually F ^ 1, a <C 1. The first layer island density is then of the order of 

iv~r-'^, (21) 

which defines the island density exponent x- In standard nucleation theory it is given by the expression 

At the moment we will however leave x unspecified to allow for the possibility of different first layer nucleation 



mechanisms. Generalizations of (t22h will be derived in Section [V 



A. Stationary nucleation regimes 

Inserting Eqs.(0|l|,|l|) into (|l|) one obtains expressions for the critical island size of the general form ||ll| 

- T^a^', (23) 

where the exponents /i and 7 are given by 

^^i^Ti' ^^wn ^'^s^'^'^ (24) 

7 = -TTT' ^ = "TTT (regime III) (25) 

7=1^^, M = (regime IV). (26) 

Consider first regime II. Consistency requires the inequalities ~ aT''' <C ix ~ a^^/*^^' and a^V^ <^ Lp ^ 
(aF)^/^ to be satisfied simultaneously, which implies 

p-(l-37)/(l-3M) <^ Q, <^ p-7/(M+l/(2»*-l))^ (^27) 

Since a ^ 1 and F ^ 1, these inequalities can be satisfied only if the exponent of F on the left hand side is smaller 
than the one on the right hand side. Inserting the expressions ( p^ ) this is found to lead to the condition 

X<^. (28) 

The analysis of regime III, based on the requirement Lx <C ic ^ -^f, yields the same condition (|2^), while in regime 
IV the relevant inequalities Ibs ^ Lc ^ Lp translate into 

X<^. (29) 

Using the expression ( |2^ ) for x, we see that the fluctuation-dominated regimes II and III become inconsistent at 
i* = i*-^ = (1 + •\/T7)/2 ~ 2.56. For larger i* regime III survives in its mean field version, while regime II disappears. 
In regime IV the transition from fiuctuation-dominated to mean-field behavior occurs at i* = 1*2 = 2; in this regime 
the critical island size is of the order of the first layer island distance, 7 = x/2, and hence the problem really reduces 
to that of first layer nucleation. 

The general form of the conditions (^,^9|) illustrates that the upper critical nucleus size, at which the mean-field 
nucleation mode sets it, is intimately linked to the expression for the first layer island density exponent. Scenarios 
diff'erent from the standard case ( p^ ) can arise if, due to e.g. metastable clusters |^ or surface reconstruction [ p^ , 
the critical nucleus on the substrate is smaller or larger than that on the island. Another example of nonstandard 
behavior is kinetically limited first layer nucleation, where Eq.(p2|) is replaced by |lj] x = 2i* /{i* + 3) and the upper 
critical nucleus sizes shift to i*i = -s/S « 1.732 and i*2 = \/TS ~ 2 w 1.6055. In addition the growth law ( p^ plays an 
important role, as will be discussed below. 
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B. Nonstationary nucleation 



For very strong step edge barriers the residence time r becomes the largest time scale in the problem, and may 
effectively be set to infinity In this situation all atoms deposited on the island remain there, and a stationary 

state where deposition and loss of adatoms compensate each other is never reached. A distinction between fluctuation- 
dominated and mean- field nucleation scenarios is nevertheless possible : In the first case nucleation occurs as soon 
as the i* + 1-th atom arrives on the island, while in the second case many more atoms have to be deposited. 

Using the growth law ( [l7| ) it is easy to show that the time required to deposit a few (i.e., i* + 1) atoms on the 
island is 

Tdcp ~ N^/yP. (30) 

At this point the island size is of the order N^^^^, corresponding to exponents 

7 = x/4, M = 0, (regime I/fi) (31) 

for the fluctuation-dominated part of regime I. To see whether nucleation is actually fluctuation dominated, the time 
scale (30) is to be compared to the encounter time (0), evaluated at the island size N~^^'^; fluctuation-dominated 



nucleation occurs if Tcnc <C Tdcp- Once more this yields the condition (|28|), which can therefore be regarded as a 
universal criterion for fluctuation-dominated second layer nucleation. For x > 2/(i* + 1) mean field theory can be 
used to compute the nucleation rate, and one arrives at the scaling exponents pl] | 

^^^^l^fr^' ^ = ° (regime I/mf). (32) 
A graphical representation of the sequence of scaling regimes as a function of i* is provided in Figure^ 



While in the stationary nucleation regimes of Section [II A the criterion (Oq) appears in a rather indirect way, in 



the nonstationary case it has a straightforward interpretation in terms of a comparison of the time scales Tdcp and 
Tenc- This is easily generalized to first layer island growth laws of the form 

L{t) - A^-i/2g„3 (33^ 

which lead to 
instead of 



IV. DIMENSIONALITY DEPENDENCE 



In statistical physics fluctuation effects are typically more prominent in low spatial dimensionalities. It is therefore 
of some interest to repeat the above considerations for growth of d-dimensional islands on a d-dimensional substrate. 
The residence time is then still given by an expression of the form (H), while the interarrival time becomes 

A^-;^- (35) 

To estimate the encounter time, it is useful to visualize the trajectories of the i* -I- 1 diffusing adatoms in a configuration 
space of dimensionality d{i* + 1), where the motion is restricted to a region of linear extent L. Nucleation occurs on 
a d-dimensional submanifold, the codimension of which is di*. We are thus dealing with a first passage problem in 
di* dimensions, where the encounter time plays the role of the mean first passage time. This implies that 



L^/v : di* <2 
L'^^* /v : di* > 2. 



Since the residence time is always at least equal to the traversal time jv^ it follows that encounter-limited nucleation 
in the sense of Tenc ^ r is possible only for dimensionalities d > 2/i*. This recovers the upper critical dimension 



dc = 2/i* for diffusion- limited reactions of i* + I particles [g^ 15 
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A. Island density exponent in d dimensions 



As a warm-up excercise we derive here a formula for the exponent x of the first layer island density valid for general 
d. The appropriate expression for the nucleation rate is ( [l6| ) , which is now to be interpreted as the rate of nucleation 
in a region of size L. The distance Id ~ N~^/'^ between first layer nuclei is determined by the condition that there 
should be one nucleation event per area ij-, during the monolayer deposition time l/F, i.e. 

u:{Id)IF ^0{l). (37) 

For di* > 2 this yields 

di* 

X = , ^ , {di* > 2), (38) 

d + 2 + 2i* ^ - ^ ^ ^ 

which agrees wit h (p2|) for d = 2 and with the expression x = j*/(3 + 2i*) derived by Kallabis, Krapivsky and Wolf 
for d — 1, i* > 2 |12||. For di* < 2 we obtain instead 

This case is physically realized only for d = = 1, where ( |39| ) reduces to the well-known value x = 1/4 p^Jl^ . 

These results are easily generalized to the growth of islands of (fractal or integer) dimensionality d' on a d- 
dimensional substrate. In that case the time of interest is not the monolayer growth time 1/F, but rather the time 
required to cover a fraction of the substrate which is of order unity. While the true coverage is related to island densitjy 
and island size through 9 ~ NL'^ , the covered fraction of the d-dimensional substrate is 9d ~ NL'^ ~ ]\[^~d./d gd/d ^ 
Setting the inverse nucleation rate equal to the time in which 6d = 0{1) is reached then yields the condition 

u:{lD)/F^f„~''', (40) 
which for di* > 2 implies the island density exponent 

di* , ^ ^ , ^ 

^ - 2(^^ + 1) + .- ^ ^''^ 

For d = 1, 2 this agrees with the general expression given in [O. 



B. Second layer nucleation in d dimensions 

For compact d-dimensional first layer islands the growth law ( p7| ) generalizes to 

L{t) - {e/Ny/'^, (42) 
and the critical island size for a nucleation rate of the form (|l^) becomes 

L^^ {QN)-^^^''+''\ (43) 
Using the expressions (OUTajlq) for the nucleation rates together with (p5|,p6|) and (Ea) it is then straightforward 



to repeat the considerations of Section [II. Rather than developing in detail the general case, we focus here on the 
transition between fluctuation-dominated and mean field nucleation regimes. 

The analysis is most transparent in the nonstationary case. One finds that the deposition time is given by ( |30| ) for 
all d, and the critical island size for second layer nucleation is 

- 7V-i/2d _ ^ _ px/2d (regime I/fl). (44) 

Using this to evaluate the encounter time (Q), it turns out that nucleation is always fluctuation-dominated for di* < 2, 
while for di* > 2 the condition for the fluctuation-dominated regime is again given by the universal relation (|28|). The 
same relation also governs the disappearance of the fluctuation-dominated regime II for general d. 
Setting the island density exponent (p8) equal to 2/{i* + 1) then yields the expression 
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i:i = ^(4-d+v/8d^ + (rf + 4)2) (45) 

for the upper critical nucleus size at which the fluctuation-dominated regimes I/fl, II and Ill/fl disappear. This is 
a decreasing function of d which approaches i*i = 1 for d — > cxd. As expected, the fluctuation-dominated regime 
becomes smaller in higher dimensions. 

To analyse the transition to mean field nucleation in regime IV we use the fact that the critical island size for 
second layer nucleation is of the order of the first layer island spacing, i.e. 7 = x/rf and Lc ^ Id- The condition for 
fluctuation-dominated nucleation therefore simply reads r/At ~ Fl]^'^ /v <C 1, or x < d/{2 + d). This clearly always 



holds for the expression (|39|), while for (38) it is true provided i* < i*2 with 

= 1 + 2/d. (46) 

Again 1*2 1 for d ^ 00, and 1*2 < i*i for all d. 

The expressions for the scaling exponents 7 and fi in the remaining regimes read 



^ x(»* + 2)-l 
d{i* + 3) ' 

X + «* 



/i = (regime I/mf) (47) 
,■* 

(regime II) (48) 



{d+l)i* + 2d' ' {d+l)i* + 2d 



M= (regime III). (49) 



i* + 2d+l i* + 2d+l 

The scaling regimes for the one-dimensional case are shown in Figure I The upper critical nucleus sizes are = 
(3 + V33)/2 « 4.372... and i*2 = 3. 

V. NUCLEATION ON THE TOP TERRACE OF A MOUND 

In multilayer growth, the suppression of interlayer transport leads to the formation of pyramidal mounds [ pT| . 
For strong step edge barriers, in the sense that /es exceeds the distance Ijj between first layer islands, the mound 
separation is set by Id and remains constant during growth |^,^ . In this growth regime the mounds are wedding- 
cake-like stacks of islands upon islands, with a characteristic up-down-asymmetry: While the valleys between mounds 
are deep crevices, at the hilltops flat terraces of a characteristic size Ltop < Id are found |jl^,|l^. Ltop is determined 
by the nucleation rate on the top terrace through the requirement that on average one nucleation event should occur 
during the growth of a monolayer, so that ||^,|iq| 

Lo[L,,^)/F = 0[1). (50) 

For a nucleation rate of the general form ( p^ this implies 

Ltop ^ r!'/'. (51) 

A more precise calculation the factor of proportionality can be found in |^ . 

As in the case of second layer nucleation, application of the expressions (^ |l5|, |l^) for the nucleation rate yields 
estimates for the top terrace size of the form 

Ltop ~ r^'a'^' (52) 
where the scaling exponents 7' and ^' depend on i* and on the scalin g regi me in question. Inspection shows that 



the equivalent of the nonstationary regime I discussed above in Section III B is trivial in the case of mound growth: 
When a <C F"^, the top terrace size becomes comparable to the lattice constant, Ltop = corresponding to pure 

statistical Poisson growth [Q. 

The fluctuation-dominated regime II with scaling exponents 

i* 

7' = = ^2 ' (regime II, mounds) (53) 
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covers the parameter range < a < T~^^, where Si = i*{2i* — l)/(2i*(i* + 1) + 2) < 1 for all i* . In marked 
contrast to second layer nucleation, here regime II exists for all i*. 

For r^*! ^ a <C T^^^^ nucleation is governed by the exponents of regime III, 

i* i* + 1 

^'^iM~3' ^' ^ IFT^ (regime III, mounds). (54) 

This regime is fluctuation-dominated for i* < 2, while for i* > 2 it splits up into fluctuation-dominated and mean 
field-like subregimes with identical scaling exponents. Both subregimes exist for all i* > 2, though the fluctuation- 
dominated part becomes very small for large i* . For a > T~^^^ one enters the regime of weak step edge barriers, where 
^ES < and the shape and evolution of the mound morphology are qualitatively different from the wedding-cake 

. A pictorial representation of the scaling regimes is given in Figure 0. 



VI. SUMMARY 

The goal of this paper has been to extend the analysis of to arbitrary critical nucleus sizes, as well as to clarify 
the relationship between the recent stochastic approaches to second layer nucleation [p|,p|pl|] and the earlier mean 
field theory of Ref. |Q. The essential new ingredient that enters the problem for z* > 1 is the encounter time Tcnc, 
which may become larger than the residence time r. The mean field approach is valid, at least as far as the values 
of scaling exponents are concerned, when r <^ Tcnc, but it fails when t > Tone- Apart from recovering the results 
of , the analysis has revealed a remarkably simple and universal criterion for fluctuation-dominated second layer 
nucleation, Eq.(p8|), which may be useful in situations where nucleation is governed by different mechanisms in the first 
and second layer. Furthermore I have verified the expectation that fluctuation-dominated nucleation should be more 
prominent in low dimensionalities, and I have shown that no upper critical nucleus size (beyond which fluctuations 
can be neglected) exists for nucleation on the top terrace of a mound. 
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FIG. 1. Scaling regimes for second layer nucleation in two dimensions. Each regime is defined by relations of the form 

p-'5fc ^ a -C r^''*', and the full lines show how the boundaries — Ina/hiF = Sk vary with the critical nucleus size i* . 
The dashed lines indicate transitions between fluctuation-dominated (fl) and mean-field like (mf) subregimcs. The explicit 
expressions for the boundaries are <5i = 1 - 3x/4 (I/fl II) ; ^2 = (x + j*)(2i* - l)/(2i*(i* -I- 1) -|- 4) (II Ill/fi); 
^3 = (3x + 2i* - 5)/{2i* - 2) (Ill/fl -» Ill/mf); ^4 = (5 -|- i* - 2x)/(6 + 2i*) (I/mf Ill/mf); ^5 = x/2 (HI IV). For 
the island density exponent the conventional expression x = «*/(«* + 2) was used. 
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FIG. 2. Scaling regimes for second layer nucleation in one dimension. The explicit expressions for the boundaries are 
5i = 1-x (I/fl -» n and I/mf-^ Ill/mf); = (x + " l)/(i* (i* + 1) + 2) (II -> Ill/fl); = 1 - 2(1 - x)/(j* - D (IH/fl 

— > Ill/mf); ^4 = X (III ^ IV). For the island density exponent x the one-dimensional expressions derived in Section IV A were 
used. 
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